Competing charge density waves and temperature-dependent nesting in 2H-TaSe2 
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Multiple charge density wave (CDW) phases in 2H-TaSe2 are investigated by high-resolution synchrotron 
x-ray diffraction. In a narrow temperature range immediately above the commensurate CDW transition, 
we observe a multi-q superstructure with coexisting commensurate and incommensurate order parameters, 
clearly distinct from the fully incommensurate state at higher temperatures. This multi-q ordered phase, 
characterized by a temperature hysteresis, is found both during warming and cooling, in contrast to previous 
reports. In the normal state, the incommensurate superstructure reflection gives way to a broad diffuse peak 
that persists nearly up to room temperature. Its position provides a direct and accurate estimate of the Fermi 
surface nesting vector, which evolves non-monotonically and approaches the commensurate position as the 
temperature is increased. This behavior agrees with our recent observations of the temperature-dependent 
Fermi surface in the same compound [Phys. Rev B 79, 125112 (2009)]. 
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The 2H polymorph of tantalum diselenide (2H-TaSe2) with 
a layered hexagonal crystal structure [1] is among the most 
studied charge density w^ave (CDW) systems. It is a quasi- 
two-dimensional metal [2] known for its incommensurate 
(ICDW) and commensurate (CCDW) density-wave phases, 
which set in below Ti^dw ^ 120 K and r^cDW ^ 90 K, respec- 
tively [3]. In its ground state, this system is a supercon- 
ductor with a rather low critical temperature, = 0.15 K 
[4]. Early neutron diffraction experiments [5] revealed in- 
commensurate Bragg peaks close to the (|00) and (|00) 
reciprocal-space positions below the second-order ICDW 
transition. Only the longitudinal incommensurability mea- 
sured along the (100) direction, was investigated, and was 
shown to decrease monotonically with decreasing tempera- 
ture, r, until it vanished abruptly at the first-order CCDW 
lock-in transition. However, x-ray diffraction measurements 
that followed shortly thereafter [6] suggested the exis- 
tence of another interjacent phase that appeared only upon 
warming at intermediate temperatures between T^cdw 
^c+ic ^ 112 K, followed by the fully incommensurate phase 
at 7 > Tc+ic- It was ascribed to a "striped-incommensurate" 
triple-q ordering, characterized by one commensurate and 
two incommensurate propagation vectors. 

In a recent angle-resolved photoelectron spectroscopy 
(ARPES) study [7-10], effects of the static and fluctuat- 
ing CDW on the temperature dependence of the 2H-TaSe2 
Fermi surface (FS) have been investigated. In particular, 
the normal-state FS has been mapped with unprecedented 
accuracy owing to its highly two-dimensional character. The 
most striking output of these experiments was the observa- 
tion of a non-monotonic T-dependence of the pseudogap 
[2, 8, 9], accompanied by a subtle, similarly non-monotonic 
variation of the FS [2, 10]. On the one hand, the per- 
sistence of these effects into the normal state, well above 
the ordering temperature, suggested that CDW fluctuations 
must be involved. On the other hand, the non-monotonic 
T-dependence could not be explained by a simple picture of 



critical fluctuations of the ICDW order that would rapidly 
decay above Ticdw. leading to a monotonic closing of the 
pseudogap. Therefore, a more intricate feedback effect of 
these fluctuations on the nesting properties of the FS could 
be suspected. This served as our motivation to take a closer 
look at the competition of commensurate and incommensu- 
rate CDW order parameters and to study the T-dependence 
of the nesting vector in the normal state. 

Being a highly accurate reciprocal-space probe of the band 
structure, ARPES can still measure the FS nesting properties 
only indirectly, e.g. via calculations of the particle-hole cor- 
relation functions based on fitting a tight-binding model to 
the measured spectra [7, 11]. The accuracy of such analysis 
is insufficient to quantify the tiny incommensurability of 
the nesting vector, which in the case of 2H-TaSe2 remains 
below 1% of the unfolded Brillouin zone (BZ) size at any 
temperature [6]. 

In the present letter, we address this challenge using 
high-resolution synchrotron x-ray diffraction. Not only does 
it provide the necessary momentum resolution for mea- 
surements of the incommensurability, but also serves as a 
more direct probe of the nesting vector due to the appear- 
ance of a thermal diffuse-scattering signal in the normal 
state [1, 12]. Here we systematically investigate its tem- 
perature dependence using x-rays with a wavelength of 
A = 0.7 A (hco = 17.7keV), characterized by a bandwidth 
of AA/A ^2x 10""^. A platelike 20-itxm-thick single crystal 
of 2H-TaSe2 was cleaved mechanically from a bigger crystal 
and mounted on a rotation stage with its c-axis parallel 
to the incident beam. The sample was then cooled down 
to 85 K with a Cry ostr earn 700 -h nitrogen cooler {Oxford 
Cryosystems) . The data were measured in transmission ge- 
ometry using a multi-purpose KUMA KM6-CH diffractometer 
and collected by successive frames with an Onyx CCD area 
detector. The orientation matrix was determined with the 
CRYSALIS software package (Oxford Diffraction, Ltd.). 

At first, we consider the temperature dependence of the 




Fig. 1 (color online). Temperature evolution of the (|00) superstructure reflection across r^cDw ^c+ic during warming (top) and 
cooling (bottom) [13]. The color represents diffracted intensity. Note the hysteretic behavior between 95 and 105 K. The momentum 
scale is given by the white hexagons, which are scaled to 1/10 of the folded BZ size. The schematic below the figure explains the origin of 
second-order incommensurate reflections in the multi-q CDW phase. Here the incommensurability is largely exaggerated for clarity. For 
an animated version of this figure with additional temperature frames, see supplementary material [13]. 



superstructure reflections across the ordered phases, mea- 
sured in the vicinity of the ( 1 0) reciprocal-lattice vector 
(see Fig. 1 or Fig. SI (a) in the online supplementary ma- 
terial [13]). In the low^- temperature CCDW phase, a com- 
mensurate superstructure reflection is observed at (|00), as 
expected. Above 87 K, it splits into five closely spaced peaks 
(Fig. 2), one of w^hich remains at the commensurate posi- 
tion, w^hereas the others represent first- and second-order 
incommensurate superstructure reflections that appear at 
(| -5, ±-^5, 0) and (| +25, ±-^5, 0) reciprocal-space 
vectors, respectively. 

The individual frames in Fig. 1 v\^ere measured at a fixed 
sample position and therefore can not accurately represent 
the full peak intensities. Therefore, we have integrated the 
diffracted intensity w^ithin several such frames, measured 
at different rocking angles, to get more reliable intensity 
distributions. The surface plot in Fig. 2 shows one such inte- 
grated pattern at T = 88 K, immediately above the CCDW 
transition, where the development of four incommensurate 
satellites can be already clearly seen, though not yet fully 
resolved. Their positions agree with those evaluated from 
Fig. 1 (b) . The origin of these multiple peaks is explained 
in Ref. 6 as a result of orthorhombic twinning of a triple- 
q superstructure with one commensurate and two incom- 
mensurate propagation vectors (C+IC), oriented at nearly 
120° to each other. The relative intensities of the peaks are 
therefore determined by the twin-domain population. The 
appearance of second-order reflections is explained in the 
schematic at the bottom of Fig. 1. 

To extract the peak positions and intensities (Fig. 3), we 
have fitted the data with a sum of two-dimensional Gaussian 
peaks. The intensity ratio between the first- and second- 
order diffraction peaks was roughly evaluated to 7 and 
was kept constant for all the fits between Tic^w and Tccpw. 
whereas the peak widths, positions and intensities were 



fitted as free parameters. Within the accuracy of our mea- 
surement, the incommensurate Bragg reflections always 
appear symmetrically with respect to the (110) direction of 
the folded BZ, so that the first-order reflections in the C+IC 
phase form an equilateral triangle with its vertex at (|00). 
Upon warming above Tc+ict ^ 107 K, the commensurate 
peak disappears, and at the same time all four incommensu- 
rate peaks collapse onto the (100) axis, losing their trans- 
verse incommensurability but preserving its longitudinal 
component, 5, without any discontinuous changes (Figs. 1 
and 3). This signifies the first-order transition to the fully 
incommensurate state that persists up to Ti^dw ^ 117 K. 
Above this temperature, broad diffuse-scattering peaks can 
be seen [13]. 

In contrast to Ref. 6, we observe the C+IC phase also upon 
cooling below T^+ic | ^ 96 K. The phase transition at T^+ic i 




Fig. 2 (color online) . Distribution of the superstructure Bragg 
intensity, measured upon warming around ( 1 0) at T = 88 K and 
integrated over several frames along the sample rocking angle. 
Note the logarithmic intensity scale. 
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Fig. 3 (color online) . (a) Intensities of the commensurate (top) and incommensurate (middle and bottom) superstructure peaks as 
functions of temperature, measured during warming and cooling. In the C+IC phase, the sum of intensities for both incommensurate 
peaks is shown. Solid lines are guides to the eyes, (b) T-dependence of the longitudinal incommensurability, 5. The inset clarifies the 
definition of 5 in the multi-q CDW phase, (c) Half width at half maximum of the diffuse peak in the normal state as a function of T. 
(d) Summary of transition temperatures between the normal state, ICDW, C+IC, and CCDW phases upon cooling and warming. 



exhibits a pronounced temperature hysteresis, in contrast 
to the lock-in transition at T^cdw that does not depend on 
the direction of the temperature change. This is best seen in 
the T-dependent intensity of the superstructure reflections 
shown in Fig. 3 (a) . The corresponding T-dependence of 
the longitudinal incommensurability and the temperature 
phase diagram are summarized in Fig. 3 (b) and (d), respec- 
tively. Note that the C+IC phase is clearly distinct from a 
simple superposition of the CCDW and ICDW phases (which 
could occur as a result of phase separation, for example) . 
Indeed, the incommensuration in the C+IC phase occurs 
along the (110) direction, whereas in the ICDW phase only 
the longitudinal component is preserved. This implies an in- 
terplay between three competing CDW vectors that possibly 
originate from the multiple local minima in the electronic 
nesting function [7, 10]. On the other hand, the fact that 
the incommensurability experiences a discontinuous jump 
at Tc+ic only in the transverse direction, whereas its longi- 
tudinal component changes continuously across the whole 
phase diagram, implies a close relationship between these 
order parameters, with a very low energy barrier separating 
the two incommensurate ground states. 

In the normal state above Tic^vv? the first-order superstruc- 
ture reflection gives way to a broad diffuse peak of nearly 
isotropic circular shape (Fig. 4), which gets weaker and 
broadens linearly with increasing temperature, as shown in 
Fig. 3 (c) . This peak can originate either from the inelastic 
scattering on the soft phonon mode [5, 12] or from the 
quasielastic scattering on the charge-density fluctuations, 
which are both peaked at the Fermi-surface nesting vec- 



tor. We have measured the positions of this diffuse peak 
in a higher BZ near the ( 1 0) wavevector to minimize the 
background intensity from the direct beam. These peak 
positions, extracted by fitting the diffraction patterns to a 
two-dimensional Gaussian profile with a linear background, 
mostly related to residual scattering from the incident beam, 
are shown in Fig. 3 (b) as circles with error bars. At Tic^vv^ 
the center of the diffuse peak coincides with the incom- 
mensurate superstructure reflection of the ICDW phase. 
At higher temperatures, however, the diffuse peak grad- 
ually shifts towards the commensurate wavevector. Hence, 
a pronounced cusp in the T-dependence of 5 is formed 
around Ticdw A similar behavior has been previously re- 
ported at the ICDW transition in SmNiCs [14]. Although 
such a non-monotonic behavior is fully consistent with the 
fluctuation-driven temperature variations of the FS observed 
in 2H-TaSe2 earlier [2, 10], it is not captured by the con- 
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Fig. 4 (color online). Temperature dependence of the diffuse- 
scattering intensity measured near the (|00) reciprocal-space 
vector. The white hexagons circumscribe the folded BZ boundary 
(i.e. are 10 times larger than in Fig. 1). For an animated version of 
this figure with additional temperature frames, see Ref. 13. 
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ventional theory of density- wave phase transitions [15] that 
considers the normal-state FS temperature-independent and 
fully determined by the static lattice potential. 

We conclude that a proper theoretical description of 
density-wave systems must take the dynamic fluctuation 
effects into account. This necessity is dictated by the ever- 
increasing accuracy of the experimental probes that can no 
longer be reconciled with a conventional (static) descrip- 
tion. As we have shown, these effects can lead to experi- 
mentally observable deformations of the FS, which, in their 
turn, exhibit a subtle feedback effect on the position of the 
nesting vector that determines the ordering instability. A 
self-consistent description is therefore required to match the 
FS nesting vector with the CDW ordering wavevector within 
a comprehensive theoretical model. 

Acknowledgements. The x-ray diffraction measurements 
reported here were performed at the Swiss-Norwegian beam- 
line BMIA of the European Synchrotron Radiation Facil- 
ity (ESRF). We are grateful to S.V Borisenko, D.V Efre- 
mov, D.V Evtushinsky A. Gabovich, V Hinkov, T. Keller, 

A. A. Kordyuk, N. Munnikes, M. Le Tacon, I. Zegkinoglou and 

B. Keimer for discussions and helpful suggestions. 



* E-mail: p.leininger@fkf.mpg.de 
^ E-mail: d.inosov@fkf.mpg.de 

[1] J. A.Wilson, E J. DiSalvo and S. Mahajan, Phys. Rev. Lett. 

32, 882 (1974); Adv. Phys. 24, 117 (1975); D. B. McWhan, 

R. M. Fleming, D. E. Moncton, E J. DiSalvo, Phys. Rev. Lett. 

45, 269 (1980). 
[2] S.V Borisenko et al, Phys. Rev Lett. 100, 196402 (2008). 
[3] The transition temperatures observed in the present study 

are consistent with those determined in the early diffraction 

experiments (Ref. 1), but are somewhat lower than those 

observed by neutron scattering (Ref. 5) . 
[4] A. M. Gabovich, A.LVoitenko, J.E Annett and M.Ausloos, 

Supercond. Sci. Technol. 14, Rl (2001). 
[5] D. E. Moncton, J. D. Axe and E J. DiSalvo, Phys. Rev Lett. 

34, 734 (1975); Phys. Rev B 16, 801 (1977). 
[6] R.M. Fleming, D.E. Moncton, J.D. Axe and G.S. Brown, 

Phys. Rev B 30, 1877 (1984). 
[7] D. S. Inosov et al, New J. Phys. 10, 125027 (2008). 
[8] A. A. Kordyuk et al, Phys. Rev B 79, 020504(R) (2009). 
[9] D.V Evtushinsky et al, Phys. Rev Lett. 100, 236402 (2008). 
[10] D. S. Inosov et al, Phys. Rev B 79, 125112 (2009). 
[11] D. S. Inosov et al, Phys. Rev B 75, 172505 (2007). 
[12] M.Holt et al, Phys. Rev Lett. 86, 3799 (2001); A.Bosak 

et al, ibid. 103, 076403 (2009); Yo Machida et al, J. Phys. 

Soc. Jpn. 73, 3064 (2004). 
[13] See Fig. SI on the next page for animated versions of Figs. 1 

and 4 including additional temperature frames. 
[14] S. Shimomura, C. Hayashi, G. Asaka, N. Wakabayashi, 

M. Mizumaki and H. Onodera, Phys. Rev Lett. 102, 076404 

(2009). 

[15] S.-K. Chan and V Heine, J. Phys. F: Metal Phys. 3, 795 (1973); 
G. Griiner, Density Waves in Solids, Addison-Wesley, 1994. 



4 



Supplementary material for the letter entitled 
"Competing charge density waves and temperature-dependent nesting in 2H-TaSe2" 
by Ph. Leininger, D. Chernyshov, A. Bosak, H. Berger, and D. S. Inosov 



(a) 



2H-TaSe, 



T=85.0K 



(b) 



2H-TaSe, 



7= 85.0 K 



warming 



CO / 



(10 0) 



(1/3 0) 



(10 0) 



{010) 



(5/3 0) 



{010) 



<<<>>>! - + 



<<<>>>! - *\* + 



Fig. SI. Animated temperature dependence of diffracted intensity in 2H-TaSe2. Use controls below the figure to start the animations, 
change the frame rate or navigate through individual frames, (a) Evolution of superstructure Bragg reflections in the vicinity of 
(|00) during warming and cooling. The momentum scale is given by the white hexagon, which is scaled to 1/10 of the folded BZ 
size, (b) Appearance of the diffuse scattering near the (|00) superstructure reflection in the normal state. Here the white hexagon 
circumscribes the folded BZ boundary, i.e. is 10 times larger than in panel (a). 



